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Gamma-rays can be produced by the interaction of a relativistic jet and the matter of the 
stellar wind in the subclass of massive X-ray binaries known as "microquasars" . The relativistic 
jet is ejected from the surroundings of the compact object and interacts with cold protons 
from the stellar wind, producing pions that then quickly decay into gamma-rays. Since the 
resulting gamma-ray emissivity depends on the target density, the detection of rapid variability 
in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be 
used to probe the clumped structure of the stellar wind. In particular, we show here that the 
relative fluctuation in gamma rays may scale with the square root of the ratio of porosity length 
to binary separation, yjh/a, implying for example a ca. 10% variation in gamma ray emission 
for a quite moderate porosity, h/a <~ 0.01. 



1 Introduction 

High-energy gamma-rays can be produced in accret- 
ing black holes with relativistic jets, as shown by the 
recent detection of a high-energy flare in the classi- 
cal X-ray binary Cygnus X-l (Albert et al. 2007). 
In this type of system the primary object is a hot, 
massive star with a strong stellar wind, and the sec- 
ondary is usually an accreting black hole. Since the 
jet propagates through the wind and the radiation 
field of the star, gamma-rays can be produced by in- 
verse Compton scattering with the radiation and/or 
pion-decays from inelastic proton-proton collisions 
with the stellar wind. If the wind has a clumpy 
structure, then jet-clump interactions can produce 
rapid flares of gamma-rays which, if detected, could 
be used to probe the size, density, and velocity of 
wind inhomogeneties. 

2 Gamma-ray emission from 
jet-clump interaction 

The basic model for hadronic gamma-ray produc- 
tion in a microquasar with an homogeneous wind 
has been developed by Romero et al. (2003). The 
reader is referred to this paper for the basic formu- 
lae. In the case of a clumped wind, the situation 
is illustrated in Figure 1. The compact object is as- 
sumed to be in a circular orbit of radius a, the clump 
crosses the jet forming an angle $ with the center of 
the star, and the viewing angle is 9. 



The spectral gamma-ray intensity (photons per 
unit of time per unit of energy-band) produced by 
the injection of relativistic protons in the clump is 

I-y(E 7 ,e)= f n(P)q^(P)d 3 P , (1) 

where V is the interaction volume, n(r') is the 

clump density, and <? T (r') is the gamma-ray emis- 
sivity, which depends on the proton spectrum and 
the interaction cross section. The spectral energy 
distribution is given by LI (Ej,0) = E 2 I 1 {E 1 ,9). 
Note that a change in the density of the target will 
be reflected in a variation of the gamma-ray luminos- 
ity. Figure 2 shows a 3D-plot of the evolution of the 
spectrum during the transit of a clump through the 
jet. For this specific calculation we have assumed 
the clump flow follows a typical (3 velocity law at 
an angle "J = 5 deg with the orbital plane, with a 
viewing angle of 45 deg. The clump density follow 
Gaussian profile with peak at ~ 10 13 cm -3 and a 
width of 0.01 R+. The primary star is assumed to 
be similar to the primary in Cygnus X-l. The frac- 
tion of the accretion power that goes to relativistic 
protons in the jet is 10~ 3 . 

From Fig. 2 we see that a gamma-ray flare with 
a power-law spectrum and a luminosity of ~ 10 35 
erg s" 1 at 1 GcV is produced by the interaction. 
The flare duration is set by the clump crossing time, 
which here is a few hundred seconds. For clumps 
interacting at higher altitudes above the black hole, 
longer timescales are possible, but the luminosity de- 
creases since the jet expands and hence the proton 
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flux decreases. 




Figure 1 : Sketch of a jet-clump interaction in a high- 
mass microquasar. 
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3 Porosity length scaling of 
gamma-ray fluctuation from 
multiple clumps 

Individual jet-clump interactions should be observ- 
able only as rare, flaring events. But if the whole 
stellar wind is clumped, then integrated along the 
beam there will be clump interactions occurring 
all the time, leading to a flickering in the light 
curve, with the relative amplitude depending on the 
clump characteristics. In particular, while the mean 
gamma-ray emission will depend on the mean num- 
ber of clumps intersected, the relative fluctuation 
should (following standard statistics) scale with the 
inverse square-root of this mean number. But, as 
we now demonstrate, this mean number itself scales 
with the same porosity length parameter that has 
been used, for example, by Owocki and Cohen (2006) 
to characterize the effect of wind clumps on absorp- 
tion of X-ray line emission. 

Let us first consider the mean gamma-ray emis- 
sion integrated along the beam. Assuming a narrow 
beam with constant total energy along its length 
coordinate z, the mean total gamma ray emission 
scales as 



Act 



n(z) dz , 



(2) 



where n(z) is the local mean wind density and a rep- 
resents an interaction cross-section for conversion of 
beam energy I\, to gamma-rays. For a steady wind 
with mass loss rate M and constant speed v ema- 
nating from a star at binary separation a from the 
z-origin at the black hole, the integral in Eq. (2) 
gives 
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(3) 



log(E 7 [BV]) 



where fi is the mean wind mass per interacting par- 
ticle (e.g., protons), and M± is the star mass loss 
rate. 

The fluctuation in this mean emission depends on 
the properties of the wind clumps. A simple model 
assumes a wind consisting entirely of clumps of a 
characteristic length I and volume filling factor /, 
for which the mean-free-path for any ray through 
the clumps is given by the porosity length h = £/f. 
For a local interval along the beam Az, the mean 
number of clumps intersected is thus A7V C = Az/h, 
whereas the mean gamma-ray production is given by 



Figure 2: 3D-plot with the light curve and the spec- 
tral energy distribution resulting from a 
jet-clump interaction. The main param- 
eters adopted are indicated in the figure. 
See text for details. 



AL ( = I(,anAz — I^anANch. 



(4) 



But by standard statistics for finite contributions 
from a discrete number AN r , the variance of this 



emission is 



I 2 a 2 n 2 Az 2 
AiVc 



I 2 a 2 n 2 hAz . 



(5) 
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The total variance is then just the integral that re- 4 PrOSpeCtS 

suits from summing these individual variances as one 
allows A 2 — * dz, 



n hdz . 



(0) 



Taking the square-root of this yields an expression 
for the relative rms fluctuation of intensity 



i 2 hdz 



f™ndz 



(7) 



As a simple example, for a wind with a constant 
velocity and constant porosity length h, the relative 
variation is just 
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TV, 
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(8) 



On the other hand, for the Owocki & Cohen 
(2006) uniform expansion model with v ~ r and 
h = h'r, we find n cx 1/r 3 and thus, 



SLy 



= Vk 



_y/j~dx/(l+xW 
J™d X /(l+X 2 ) 3 / 2 



y/2h'/3. (9) 



Typically then, if h ~ 0.03a, 
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0.1. 



This implies an expected flickering at the level of 
10% for a wind with such porosity parameters. The 
variability timescale will depend on the wind speed, 
the size of the clumps, and the width of the jet, with 
a typical value of ~ an hour. 



Fluctuations of 10% in a source with a luminosity 
of 10 34-35 erg s _1 on timescales of ~ 10 3 s could 
be detectable with GLAST and CTA if the source is 
located at a few kpc. This means that gamma-ray 
astronomy can be used to probe the structure of stel- 
lar winds through dedicated observations of micro- 
quasars with massive donor stars. If additional infor- 
mation on the wind can be obtained through X-ray 
observations (Owocki & Cohen 2006) and the source 
can be detected with neutrino km-telescopes, we 
can also gain valuable insights into the relativistic 
proton content of the jets. 
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